2D Stimulated Resonance Raman Spectroscopy of Molecules with 
Broadband X-ray Pulses 

Jason D. Biggs, 1 Yu Zhang, 1 Daniel Healion, 1 and Shaul Mukamel 1 '^ 
Department of Chemistry, University of California, Irvine 
Irvine, CA 92697-2025, USA. 

(Dated: 3 March 2013) 

Expressions for the two-dimensional Stimulated x-ray Raman Spectroscopy (2D-SXRS) signal obtained us- 
ing attosecond x-ray pulses are derived. The ID- and 2D-SXRS signals are calculated for trans-N-methyl 
acetamide (NMA) with broad bandwidth (FWHM ~14.2eV, 181 as) pulses tuned to the oxygen and nitrogen 
K-edges. Crosspeaks in 2D signals reveal electronic Franck-Condon overlaps between valence orbitals and 
relaxed orbitals in the presence of the core hole. 



I. INTRODUCTION 

Driven by the development of highly intense attosec- 
ond sources of x-ray radiation^ there is much inter- 
est in mapping nonlinear spectroscopic techniques de- 
veloped for the NMR, IR and visible regimes to x-ray 
frequencies PS These techniques use sequences of laser 
pulses with well-defined inter-pulse delays to probe the 
correlation between different nuclear spin, vibrational, 
or electronic molecular eigenstates. One-dimensional 
stimulated Raman s-ray spectroscopy (1D-SXRS) has 
been proposed to prepare and probe valence electronic 
wavepackets during a delay between two pulses.™ We 
present the two dimensional stimulated x-ray Raman 
spectroscopy technique (2D-SXRS), an extension of 1D- 
SXRS which may be used to probe molecular valence 
electronic excitations in greater detail. First we briefly 
survey the optical domain, time-resolved, stimulated Ra- 
man experiments which inspire these new x-ray tech- 
niques. We then will discuss the dynamics probed by 
the time and frequency domain x-ray spectroscopies, in- 
cluding 2D-SXRS. 

Resonances are observed in optical or ultraviolet Ra- 
man spectroscopy when the difference frequency uii — u>2 
of the incident (wi) and scattered (012) field modes match 
a vibrational transition frequency!^ This inelastic scat- 
tering process can be either spontaneous or stimulated 
by a second interaction with the applied field!^] Fem- 
tosecond lasers with bandwidths comparable to molecu- 
lar vibrational frequencies (on the order of hundreds of 
wavenumbers), have made it possible to excite vibrations 
impulsively!^^ The dynamics of a slow vibrational sys- 
tem are probed through a perturbation that depends 
on its coupling to a fast electronic systemP ° l 14 l 15 l Time- 
domain vibrational Raman spectroscopy with picosecond 
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delays between pulses was developed in the seventies 9 
and used to probe vibrational dephasing. The impulsive 
pump-probe technique^ ! 13 ! 16 ! is one dimensional (ID) 
since a single delay time is controlled. Selecting a spe- 
cific intermediate electronic state by tuning the central 
frequencies of the exciting pulses determines a unique ef- 
fective nuclear perturbation during the Raman process. 
Multidimensional time-domain stimulated Raman spec- 
troscopies extend the stimulated Raman technique to 
multiple time evolution periods. Fifth-order off-resonant 
stimulated Raman spectroscopy was first proposed to ex- 
amine the vibrational structure in liquids.^ 

These nonlinear techniques were developed in paral- 
lel with technological progress in the coherence and in- 
tensity of ultrafast lasers. A new generation of very in- 
tense ultrafast x-ray sources allow weak nonlinear x-ray 
transitions to be observed for the first time despite the 
low cross-section and short lifetimes of core-excitations. 
X-ray free electron lasers (XFELs)P5i convert the kinetic 
energy of a beam of accelerated electrons to an electro- 
magnetic field by passing the beam through a magnetic 
grating.ES The Linac Coherent Light Source (LCLS) at 
Stanford 20 ! 21 ! j s a nex t generation XFEL capable of gen- 
erating extremely intense (~ 10 13 photons) pulses with 
x-ray central frequencies (560 — 2660 eV), 1 and has al- 
ready been used to probe exotic forms of plasma 22 and 
to find the structure of biological samples from a suc- 
cession of single protein x-ray diffraction snapshotsJ^Hil 
The source is intense enough to create hollow atoms by 
depopulating both core orbitals for a resonant transition 
within the Auger decay time of the core hole in atomic^ 
and molecular^ systems. One drawback of XFEL radi- 
ation is that the coherence time is short (< 5 fs),^ a 
consequence of quantum statistics in the electron bunch 
used to generate the pulse. Despite this drawback, the 
high intensity of the LCLS XFEL source make it the most 
likely candidate for exciting nonlinear x-ray transitions. 
The technology needed to create pulses with desired char- 
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acteristics is rapidly improving; an XFEL pumped x-ray 
laser was recently demonstrated, capable of generating 
intense, temporally coherent x-ray pulses.®! Similar ad- 
vances promise to significantly extend the capability of 
the LCLS light source. 

Pump-probe and stimulated Raman techniques were 
among the first nonlinear experiments performed in 
the visible regime Ehfil2flJ This paper focuses on x-ray 
analogues of multidimensional stimulated Raman tech- 
niques, using soft or hard x-ray pulses with bandwidths 
greater than 10 eV to create wave packets of electronic ex- 
citations. The x-ray interaction involves single particle, 
field driven, transitions between core and valence elec- 
tronic orbitals and the many-body valence response to 
a transiently created core-hole. 30 Extensive theoretical 31 
and experimental 3 !] work was applied toward using this 
unique source as a probe in time-dependent experiments. 
The low cross-section of the x-ray matter interaction 
make IR and optical pump/x-ray probe experiments 
more accessible than the all- x-ray pump-probe technique. 
IR pump/x-ray probe experiments were proposed and 
simulated for diatomic molecules JHHH1 Time-resolved op- 
tical pump/x-ray probe experiments 37 have been per- 
formed at the LCLS which highlight the experimental 
difficulties in synchronizing pump and probe pulses. 32 
All-x-ray stimulated Raman experiments in which the 
pump and probe pulses both have x-ray frequencies are a 
natural extension of these techniques. Frequency-domain 
resonant inelastic x-ray scattering (RIXS) is a well estab- 
lished technique that probes the single-particle occupied 
and unoccupied density of states around the resonant 
core by measuring the energy-resolved x-ray radiation 
spontaneously emitted by a core- excited state l^Mi An 
x-ray photon is absorbed, exciting a core electron to an 
empty (virtual) orbital, and an electron drops into the 
core-hole, emitting a photon. Peaks in the Fourier trans- 
form of the time-domain ID- SXRS^ signal represent ex- 
citations of valence electronic states generated through 
a core-excited intermediate. Ultrashort x-ray pulses are 
difficult to manipulate 4 ^ and the cross-section for core- 
excitation is lower than for excitations at lower frequen- 
cies. The feasibility of stimulated and coherent Raman 
experiments using the LCLS was explored, 42 and schemes 
for tailoring the electron bunches to generate pairs of at- 
tosecond pulses suited for stimulated x-ray Raman ex- 
periments have been suggest 

Simulations® and experiments^ suggest that electron 
correlation can drive charge migration in ionized elec- 
tronic systems. Thanks to the elemental specificity of 
resonant x-ray excitation, it can be used as a perturba- 
tion that is local in time and space, making nonlinear x- 



ray spectroscopies attractive for measuring ultrafast dy- 
namics with high spatial resolution. Nonlinear effects in 
the propagation of XUV pulses,^ and atomic attosecond 
electronic dynamics, 47 were investigated in noble gases, 
and show that nonlinear effects contribute to measure- 
ments using XFEL sources. The 1D-SXRS signal can 
be expressed as the time-dependent overlap of an im- 
pulsively excited valence electron doorway wavepacket 
with a stationary window created by the probe. SXRS 
was proposed to investigate the electronic properties of 
molecules,^ and applied to conjugated 7r-bonded or- 
ganic molecules 4 and excitonic systems^ using a tight- 
binding model Hamiltonian. Several extensions and re- 
finements of 1D-SXRS were proposed, including the use 
of attosecond pulses to prepare entangled particle-hole 
states, 4 8J frequency domain coherent anti-Stokes Raman 
spectroscopy with wide- and narrow-band pulses 50 } and 
a many-body Green's function-based method suitable for 
calculating the SXRS of larger systems.^ X-ray spectra 
of molecules with cores separated from each other at some 
distance were discussed™ 

This paper extends the formalism previously applied 
to 1D-SXRS to two dimensions (2D-SXRS) by adding one 
more pulse. Expressions for the 1D-SXRS and 2D-SXRS 
signals are presented in sections [IT] and |III| respectively. 
Simulations of the UV, ID and 2D-SXRS signals for the 
nitrogen and oxygen K edges in trans-NMA are reported 
in Sec. IV and their significance is discussed in Sec. [V] 



II. ULTRAFAST ONE DIMENSIONAL STIMULATED 
RAMAN: 1D-SXRS 

1D-SXRS is the simplest time-domain, x-ray Raman 
technique. The electric field is represented as 



E(r, t) =ei£i(t) exp(iki • r — iwit) + e2£2{t — t) 
x exp(tk2 • r — iu2(t — r)) + ex.. 



(1) 



Here kj, Wj, ej and £j are the wavevector, carrier fre- 
quency, polarization vector, and complex envelope of the 
jth pulse, respectively. The system is excited by a pump 
pulse and the transmission of the probe pulse is recorded 
after a delay t (see Fig. [I) ] 4 ' 52 ' 53 ' 

Each pulse interacts with the system twice via a stim- 
ulated Raman process. A core-electron is excited into an 
unoccupied orbital during each pulse; the excited system 
then evolves for a short period before a second interac- 
tion with the same pulse stimulates a valence electron 
to destroy the core-hole, emitting an x-ray photon in the 
process. The system is left in a coherent superposition of 
valence excited states after interaction with the first x-ray 
pulse which interferes with the amplitude of the Raman 
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FIG. 1. Loop diagrams, pulse sequence, and core and valence energy levels for the 1D-SXRS signal. 



process with the second pulse, leading to a change in its 
transmitted intensity.^] The 1D-SXRS signal is defined 
as the change in the transmission of the second pulse 
with or without the pump pulse, as a function of the de- 
lay between the two pulses. We neglect nuclear motions, 
further work will be needed to examine their effect on 
SXRS. 

We consider stimulated resonance Raman excitation 
at soft x-ray wavelengths, which are much longer than 
the spatial extent of the resonant core orbitals. A phase 
dependent on the distance between resonant cores will 
not enter into the response, as each field interacts twice, 
canceling the phase factor of exp(±ikj • R n ). Taking the 
delay between pulses (r) to be long compared with the 
lifetime for core-excited states (< 10 fs for nitrogen and 
oxygen), 54 the signal will be dominated by the ground- 
state bleach (GSB) contributions to the time-integrated 



pump-probe signal (defined as the pump-induced change 
in the absorption of the probe). 

The closed-time path loop diagrams for the two in- 
terfering contributions to the signal are shown in Fig [T] 
These diagrams are similar to the double-sided ladder 
diagrams traditionally used to describe nonlinear optical 
spectroscopy, 5 with a few key differences. We read the 
loop diagram by starting with the ground state on the 
lower left corner and moving upwards, describing forward 
evolution (propagation of the ket). Interactions with the 
field are represented by arrows facing into (absorption of 
a photon) and out of (emission of a photon) the diagram. 
At the top of the loop we switch to backward time evo- 
lution (propagation of the bra), and finally arrive back 
at the ground state. Thus we are able to work in Hilbert 
space, rather than the higher dimensional Liouville space 
required with the density matrix (for diagram rules see 
Refs. HHandHU. From Fig [I] we get 
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(2) 



+ c.c. 



Here \g) is the electronic ground state, |e) is a state with where 

a core hole, and \g') is a valence excited state (see Fig. .oo r r 2 

[IJ. Eq. [2]can be recast as a j;g'g" — ~ ^^( e j * ^g'e)( e j * Veg") I dr 2 I dr\ 

e J —CO J —CO 

Ssxhs(t) = 0? [(a 2 (r)a 1 (0)) - (a\ (0)a 2 (r))] (3) x S * ( T2 ) £j ( n ) x eX p (iti eg ,T 2 - iA^„7i) 

(4) 
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is the effective polarizability weighted by the two-photon 
spectral density of the j th ultrashort pulse. Here A : ev = 
Wj — w e v + iF e is the detuning for the v — » e transition, 
and T e is the inverse of the core-hole lifetime. In Eq. [4] 
we have included the direction cosines between the pulse 
polarization vector and the transition dipoles, as well as 
the core-excited lifetime Y e . Eq. [4] can be recast in the 
frequency domain as 



a r,g'g" ~ E 



[ej ■ V g , e )(ej ■ V eg n) 



2tt 



x / dw 2 



(5) 



"2 + &lg> 



The first term in Eq. [3] (diagram I from Fig. [TJ can be 
viewed as a valence wavepacket cti \tpo), created by pulse 
1, which propagates forward in time r and overlaps with 
a window wavepacket (^o| a 2 created by pulse 2. The 
second term (diagram II) can be viewed as a window 
wavepacket created by pulse 2 c*2 | Vo) propagating back- 
wards in time — r to overlap with the doorway (^o| «i 
created by pulse 1. 

Hereafter we assume linearly polarized Gaussian 
pulses, with the spectral envelope function 



(6) 



where tTj is the temporal pulse width (equal to intensity 
FWHM divided by 2v^n2). The absolute phase <pj of 
the j th pulse does not affect the signals considered here. 
Inserting Eq. [6] in Eq. [5] the effective polarizability 



becomes 



a j;g'g" ~ ~ ^ 



Y^(e d ■ V g , e )( ej ■ V efl «)£i (A J afl ,)£i(Ai fl „) 

e 

/ \% + ei&{-<T j (Ai gl + Ai g „)/2)] . 

(7) 



where 



erfi(a;) 



2i 



V* Jo 



die' 



(8) 



is the imaginary error function (see Appendix [A] for de- 
tails). The treatment so far (Eqs. [4j [5] and [7) has 
assumed the molecule is oriented in the lab frame. Of 
course, under typical experimental conditions in the gas 
or condensed phases, the material will be an ensemble 
of randomly oriented molecules. All signals presented 
below are orientationally averaged, assuming all incom- 
ing fields are polarized parallel to each other, using the 
framework of Ref. [57) The necessary tensor expressions 
for the 2D-SXRS signal are given in Appendix [B| 



The form for the effective polarizability given in Eq. [5] 
is general, and can easily be extended to cover the possi- 
bility that the upward and downward transitions are fa- 
cilitated by different pulses whose arrival times coincide. 
The theoretical formalism presented is also well suited 
explore the use of pulse shaping^"^ to optimize the Ra- 
man signals and highlight desired features. The experi- 
mental difficulties for x-ray pulse-shaping are daunting, 
but technological progress in this area is proceeding very 
rapidly. 

The operator aj is complex and symmetric and there- 
fore non-Hermitian. 62 We define its Hermitian and anti- 
Hermitian components as 



af =(ay + a})/2 
af H =(a,-a])/2 



(9) 



The matrix elements of these operators are easily found 
from Eq. [7j 



a j;g'g" ~ 2 



e 

xerfi(-^(A^+Ai g „)/2) 



The 1D-SXRS from Eq. [3] may be rewritten as 

Ssxfls(r) = 2X{at H (T)a 1 (0)). 



(10) 



(11) 



This signal depends on only on the anti-Hermitian part of 
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the polarizability for the probe pulse, which ,as Eq 
shows, decreases more rapidly with off-resonant detun- 
ing. We may ignore vibrational contributions to the elas- 
tic component of the signal which depends on oii- gg . The 
short core-excited state lifetime precludes a vibrational 
phase evolving on this potential energy surface; the x- 
ray polarizability is diagonal in the vibrational subspace 
within the Condon approximation. Vibrational progres- 
sions may appear for long delays as an additional fine 
structure to the ID and 2D signals. These are not in- 
cluded in the present simulations. 



III. TWO DIMENSIONAL STIMULATED RAMAN: 
2D-SXRS 

In 2D-SXRS the transmitted intensity of the third 
pulse depends on two pulse delays. The 2D signal is 
related to x^ h \ whereas the 1D-SXRS pump-probe sig- 
nal is associated with x( 3 ).E^The extra time-delay period 
allows more complicated valence electronic dynamics to 
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be prepared and probed. Again we assume that the delay 
times ti and t 2 are long compared to the core-hole life- 
time, any core-excited populations created by the pulses 
may be safely neglected. The 2D-Raman signal is repre- 
sented by the four diagrams in Fig. [2] Using a single- 
particle picture, the first pulse creates an electron-hole 
pair in either the ket (diagrams i and iv) or the bra (di- 
agrams ii and Hi), The second pulse can either create 
another electron-hole pair (as in diagrams i and Hi), or 
change the electron-hole pair created by the first (dia- 
grams ii and iv). In either case the pairs created by the 
first and the second pulses must share either a common 
hole or electron in order to survive the trace. These di- 
agrams represent different sequences of forward and/or 
backward time evolution of valence wavepackets. These 
generalize Fig. [I] which has a single evolution period. 
The expression for the Raman signal can be read directly 
from the diagrams in Fig. [2] and is given by 

S2D-sxR S {ti,t 2 ) = 9 [Si + S tt + S ttt + Si„] , (12) 
where 

S* = -<aS(*iW*a+tOai(0)) 
S«H«I(0)o4(*i)aa(*2 + ii)) (13) 
S«i = -<«*! (0)a 3 (ta+ti)a a (ti)) 
Si„ = (a 3 (ia+*iW*Oai(0)) 



This signal may also be recast as 
®2D-sxR S {ti,t 2 )=2S[{af H {t 1 +t 2 )a 2 {t 1 )a 1 {0)) 

The first term in Eq. [l4j Sa + is the time-dependent 
overlap, within a resonance window determined by the 
probe pulse, between the ground state and a valence 
wavepacket that has interacted with both pump pulses. 
The second term, Sj + Sm, is the overlap between differ- 
ent Raman wavepackets. 

IV. SIMULATIONS FOR THE OXYGEN AND 
NITROGEN K-EDGES OF NMA 

A. Electronic Structure Calculations 

The trans-NMA core-excited states are modelled as 
single core- hole/ virtual orbital electron pairs in the 
Static Exchange (STEX) molecular orbital (MO) ba- 
sis using the orbital approximation. Valence excited 
states are treated at the Configuration Interaction Sin- 
gles (CIS) level of theory. Both in STEX and CIS 
Hartree-Fock orbitals are employed. In order to keep 
the level of theory consistent we have calculated the 
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excited states using CIS, which is known to overesti- 
mate the adiabatic excitation energies.^ More accurate 
TDDFT methods to calculate the core excitation energies 
are under development.^ The geometry of irans-NMA 
was optimized using the quantum chemistry package 
Gaussian09^ at the B3LYP/6-311G** level. The STEX 
mo del EsHU] m which a core-ionized Fock matrix is con- 
structed which includes the core-hole self-consistently, 
was implemented in a modified version of the quantum 
chemistry package PSI3, 69 with an implementation de- 
scribed in Appendix C of Ref. 53j All transition fre- 
quencies and dipole moments for core-excited states were 
calculated using the orbital approximation at the HF/6- 
311G** level. Fig. [3] shows the largest amplitude Natu- 
ral Transition OrbitalsP2l(NTO) between the ground and 
excited states which contribute to Figs. |9~||T3"1 obtained 
by a singular value decomposition of the CIS transition 
densities. 

Using Fermi creation and annihilation operators for 
the valence (c) and core (d) molecular orbitals, the core- 
excited states are given by 

|e) = \jn) = c]d n \e Q ) (15) 

where |e ) is the lowest energy core excited state. The 
valence excited states are 

l5') = E C7 -l«) = E C7 -4^3)- (16) 

ai ai 

1 

where C ai are the CI coefficients for state \g'). 

B. Valence Excitations and UV Absorption 

Earlier simulations on trans-NMA presented in Ref. 
[531 approximated the valence excited states as single 
electron-hole pairs in the HF MO basis. The lowest 
energy valence-excited state then involved a transition 
from the highest occupied orbital to the lowest virtual 
orbital, a HOMO to LUMO transition. The excitation 
energy was approximated as the difference in MO en- 
ergies, which gave 14.35 eV. Here we use a higher level 
treatment of valence excited states obtained by diagonal- 
izing the singly-excited block of the Hamiltonian. The 
lowest-energy excitation is now found to be 6.92 eV (see 
Table [ijfor a list of the CIS energies), which is consistent 
with experimental vacuum-UV absorption results from 
trans-NMA:2B The valence excitations and the Raman 
peaks presented here (between 7 and 18 eV) are very 
different from those in Ref. [53] 

The calculated UV absorption spectrum is displayed 
in Fig. [4] The molecule is taken to be initially in the 
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ground state and the 8{w — cv g i g ) in Fermi's golden rule 
was replaced by a Lorentzian, with an effective lifetime 
Y g i, to account for the linewidth of valence electronic 
transitions 

An accurate estimate of Y g i would require potential en- 
ergy surface calculations and characterization of any con- 
ical intersections or dissociative states which could cause 
decay of the excited state population in \g'). We set 
Ygi = 0.05 eV for all valence excitations in all simula- 
tions for display purposes. 



C. Frequency Domain XANES and RIXS Spectra 

In Fig. [5] we present the x-ray Absorption Near-Edge 
Sstructure (XANES) spectra calculated for an isotropic 
trans-NMA sample, excited at the nitrogen and oxygen 
K-edges. The XANES and the UV absorption spectra 
were calculated using a Lorentzian lineshape to account 
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FIG. 3. Natural transition orbitals of the dominant excitations in Fig. |llj|13| 
for core lifetime broadening and pure dephasing 

g pjS&W <18 » 



As in Ref. 53, we set the first XANES transition for ni- 
trogen to 401.7 eV and to 532.0 eV for oxygen, to match 
experiment. 72 Each XANES peak represents an excita- 
tion from the core-orbital to a bound virtual orbital. Or- 
bitals with energies above the ionization potential are 



coupled to a continuum of photoelectron states. The 
linewidth for nitrogen was taken as: = 0.4 eV for 

energies up to 408 eV, and ramped up to 1.5 eV at 415 
eV and held constant for higher energies. A similar form 
was followed for oxygen: Y e o = 0.4 eV for energies up 
to 537 eV, and ramped up to 1.5 eV at 544 eV and held 
constant thereafter. These phenomenological linewidths, 
which were used in earlier studies to match experimental 
XANES spectra of small nitrogen and oxygen containing 
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FIG. 5. Simulated XANES from trans-NMA at the nitrogen 
(top) and oxygen (bottom) K-edge. The stick spectra (black 
lines) have been convoluted with a lineshape function (see 
text) to give the spectra in red. Shown in blue are the power 
spectra for the Gaussian pulses used in the time-domain ex- 
periments described here. 



organic molecules J 5 2 1 73 1 can reflect a variety of broadening 
mechanisms, including vibrational motion of the core- 
ionized molecule, autoionization of the excited electron 
and direct photoionization of the core into a photoelec- 
tron state. An exact treatment of these effects would re- 
quire a Wigner-Weisskopf treatment of the decay of the 
bound excited electron into a continuum of photoelec- 
tron states, and potential energy surface calculations of 
these core-excited states. We use the atomic core- lifetime 



given in Section [TT| in the calculation of a. 

To set the stage we first calculate the traditional 
frequency-domain spontaneous Raman (RIXS) spectra. 
These should serve as a reference for comparing with 
the 1D-SXRS signals. The optical resonant Raman tech- 
nique probes which vibrational modes of a molecule are 
perturbed by a given electronic excitation. Only those 
modes whose potential is different in the ground and ex- 
cited electronic states have Franck-Condon activity. In 
RIXS a core electron is excited into an empty virtual 
orbital and then de-excited. Inelastic losses represent- 
ing valence excitations can be created directly, when an 
electron other than the excited core electron drops into 
the core-hole, or indirectly when the Coulomb poten- 
tial of the transiently-created core-hole creates valence 
excitations) 74 ! 75 ! The energy satellite spectra for the in- 
direct process represent valence shake-up states, and are 
dependent on the excitation frequency. 

In RIXS a monochromatic x-ray beam (w{) impinges 
on the molecule, and the scattered radiation (w 2 ) is fre- 
quency resolved, with peaks recorded vs. u>i — w 2 = <*> g i g . 
Here we take the incident beam w\ to be polarized, and 
the scattered light w 2 is also sent through a polarization 
filter before detection (detected parallel to the excitation 
polarization). The Kramers-Heisenberg formula with a 
Lorentzian, linewidth Y g i =0.05 eV, gives 



^RIXS 



(wi,w 2 ) 



E 

s' 



l 9'9 



("i)l 2 rv 



(Wj - w 2 



w g , g y + r 2 . 



where 



E 



(e/ -/VeXe/ -Meg) 



u eg + tr e 



(19) 



(20) 



is the frequency-dependent resonant polarizability. Eq. 
[20] applies to an oriented molecule in the lab frame. To 
account for the isotropic distribution of molecules, we 
orientationally averaged Eqs. [l9| and |20] assuming paral- 
lel excitation and detection polarization.^ 

The calculated RIXS spectra with the incident light 
tuned to the nitrogen and oxygen K-edge are displayed 
in Figs. |6]and|7] respectively. The RIXS intensity falls off 
as the laser frequency is detuned from the strong core- 
edge transition, as expected from Eq. |20j The peak 
pattern changes with Wi, since the intensity also depends 
on the transition dipole between the valence-excited state 
in question and the core-excited state in resonance with 
the laser frequency. 

When Wi = 401.7 eV, the transition frequency for the 
first Nls core-excited state, the valence-excited states S 2 
and 55 are strongly Raman active. Tuning Wi to be 
resonant with the second core-excited state causes these 
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FIG. 6. Calculated RIXS signal at the nitrogen K-edge from 
trans-NMA. 

peaks to be replaced by a strong S3 peak. Changing the 
excitation frequency to 406.5 eV (the ninth core-excited 
state) we see the absence of S 2 , S 3 , and S 5 peaks, and the 
dominant peak is now from the S lg state. The oxygen 
RIXS spectrum, shown in Fig. [7] shows a similar pattern. 
The signal is strong at the core edge, Wi = 532.0 eV, and 
is dominated by the S 3 , S 15 , and S lg peaks. Increasing 
the excitation frequency to 536.6 eV gives strong peaks 
for Sn and S27. This atomic selectivity makes RIXS a 
powerful tool for probing the derealization of valence 
excited states. 17 - 7 - 17 ^ 
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FIG. 7. Calculated RIXS signal from trans-NMA at the oxy- 
gen K-edge. 



While the valence-excited states that contribute to the 
UV absorption spectrum and the RIXS and SXRS spec- 
tra are the same, their intensities are radically different, 
as can be seen by comparing Pig. [4] with Pigs. [6] and 
[7] This is because the absorption is determined by the 
dipole moment between the HF ground state and the 
valence-excited state, while the Raman spectra depend 
on the polarizability resonant with a given core-excited 
state. 

D. Time-domain Stimulated X-Ray Raman Spectroscopy: 
1D-SXRS 

In both 1D-SXRS and RIXS 79 , information on valence 
excited states is obtained through transient excitation of 
core excited states. 1D-SXRS has two advantages over its 
frequency domain counterpart. First, because the inter- 
action times are constrained by the field pulse envelopes, 
it should be possible to probe electronic dynamics di- 
rectly by preceding the two Raman pulses with an ini- 
tiation pulse. This pulse could prepare the molecule in 
some nonstationary state which may then be probed by 
an SXRS experiment, in much the same way as is cur- 
rently done in vibrational spectroscopy!^^ We do not 
pursue this course here, rather we look at another advan- 
tage inherent in the time-domain experiment. We can 
use a two color scheme, where the two pulses are tuned 
to be resonant with different core transitions, providing 
an additional experimental knob to turn. The differences 
between SXRS spectra with different pulse configurations 
could help shed light on the nature of the valence excited 
states of a molecule. The two diagrams in Fig. [I] recast 
the SXRS spectrum as a two-slit experiment where va- 
lence wave packets created on either the ket or the bra by 
the first pulse are probed by another wave packet created 
in the ket by the second pulse. 

The calculated 1D-SXRS signal is displayed as the 
Fourier transform of Eq. [3] 



5>SXHS 



(0) 



p OO 

: J drEsx rs( t ) exp [iQr] 

D? {a2;gg'Otl-g'g + (<* l;gg> )*a 2 ;g'g) g' - tfl) + S {0-1;ggi 0-l;g> g ~ {<X\;gg> )* 0-2\g' g) Ug'g 

(rf, -2tr fl ,n-n2 + w 2g 



Unlike RIXS (Eq. [l9]), §sxrs(^) results from the in- 
terference of the two pathways in Fig. [I] and may not 
be written as an amplitude squared. The signal is col- 
lected in the time domain, and Fourier transformed nu- 
merically. It is important that the interpulse delay be 



longer than the core-hole lifetime (> 10 fs). In this case, 
any interaction with the pulses that results in a core- 
excited population will decay to an ionized state through 
an Auger process, with transitions that are detuned from 
the resonant excitation with the probing pulse. If the 
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core-hole lifetimes were long compared to the delay, the 
signal would be distorted by an interference between 
the wavepacket produced on the core-excited states and 
the ground valence electronic wavepacket. 52 We also ne- 
glect a r-independent contribution to the signal, result- 
ing from elastic scattering with g' = g. This can be 
removed experimentally. 
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valence-excited state number 

FIG. 8. Hermitian and anti-Hermitian parts of the effective 
isotropic polarizabilities (Eq. [7j for the two pulses used in our 
simulations corresponding to the nitrogen and oxygen K-edge 
excitations, plotted using an arcsinh nonlinear scale (shown 
on the right). The Hermitian part is purely real, while the 
anti-Hermitian part is purely imaginary. The axes are labeled 
by the state numbers, for the ground state, 1 for Si, etc. 
State assignments can be found in Table [I] 



The simulated 1D-SXRS signals are shown in Pig. [9] 
Parallel field polarizations are used, and the signal is av- 
eraged over an isotropic distribution of molecules.^ The 
4 rows represent the four possible pulse configurations, 
e.g. row 3 is labeled ON signifying that the first pulse is 
tuned to the oxygen and the second pulse to the nitro- 
gen K-edge. The three columns show the real, imaginary, 



and modulus of Eq. 21 



The one-color, time-domain signals in the top two rows 
of Fig. [9] resemble their frequency-domain counterparts 
of Figs. [6] and [7] In the N-edge RIXS spectrum, the 
S2 and S3 peaks are strong when the excitation fre- 
quency is set to 401.7 eV and 402.5 eV, respectively. 
Both peaks show up prominently in the NN 1D-SXRS 
spectrum. Since the excitation must be broadband to co- 
herently excite valence excited states, 1D-SXRS does not 
select a given e state. An interesting possibility will be to 
use a combination of narrow band excitation combined 
with broadband stimulated de-excitation, as is commonly 
done in optical Raman spectroscopy,^ to regain this se- 
lectivity over core-excited states. 

Two-color 1D-SXRS signals have no RIXS analogue. 
In the third and fourth rows of Fig. [9] we show the ON 
and NO spectra. The moduli of these two signals are 
virtually identical but differ in their real and imaginary 
parts. The S 2 peak is mostly absent from the OO signal, 
very strong in the NN signal, and moderate in the ON 
and NO signals. However, the two-color signals are not 
simply averages of the one-color signals from the top two 
rows. For instance, the S10 peak that shows up promi- 
nently in the NN spectrum is missing completely from the 
two-color spectra. Thus, two-color 1D-SXRS provides a 
different window into the interaction between core-hole 
excitation and valence-excited states. 



The Hermitian and anti-Hermitian parts of the res- 
onant polarizability matrix a averaged over the pulse 
bandwidth are shown in Fig. [8] We assume pulses 
with Gaussian envelopes, durations of <Tj = 77 as ( 
1 /<7j ss 8.5 eV), and center frequencies set to the N and O 
K-edge transitions (401.7 eV and 532.0 eV, respectively). 
The axes are labeled by the valence excited state num- 
ber, with referring to the ground state. We display the 
isotropic polarizability dj, obtained by replacing the di- 
rection cosine term (e.j ■ V g i e ){ej ■ V eg ") with V g i e ■ V eg n 
in Eq. |4j Otherwise, it would be necessary to adopt 
a given molecular lab-frame orientation to visualize aj. 
This ID SXRS spectrum probes matrix elements of the 
effective polarizability in Eq. [7]between the ground state 
and the set of singly excited state, aj- gg i for j = N, O. 
The 1D-SXRS signal is a function of the top row of this 
symmetric matrix. 



Two-dimensional stimulated Raman: 2D-SXRS 



The 2D-SXRS signal in Eq. 13 is the product of o-i <g ' g 
and ak- gg " , from the top row of the polarizability matrix 
with a third term, a; ; g' fl » from its interior (see Fig. |8j. 



The 1D-SXRS signal (Eq. (21| only depends on the first 
row and column of the polarizabilities; 2D-SXRS thus 
carries information about correlations between dynamics 
in the two time periods not available in the ID signals. 

The 2D signals will be displayed in the frequency do- 
main 



^2D-SXRS(&1, ^2) = / dt 1 

Jo 



\ l df 2 e initl+in2 ' 2 8 2D -sx.Rs(*i,*2), 

(22) 



11 



Re 



Im 



Abs 



• i— i 

a 



• i— i 

Pi 

CD 



0.5 
ILK 
-0.5 



i.l.IKi 
-0.25 
-0.50 
-0.75 
-1.00 
-1.25 

-1.50 

I. 00 

0.75 
O.SO 
0.25 
0.00 
-0 25 
-0.50 

i on 

0.75 
0.50 
0.25 
0.00 
-0.25 

II. 50 



NN 






/r"V 






s 

ft fto 


19 

llftl 


OO 






r 


■Sis 

c ^ 3 */ 

i JV " 


Sl9 
jftl 


ON 










ft 1 


ft 

Sis 

I 


Sis 


NO 


^\ 








_U 


ft B 



10 12 14 16 6 



10 12 14 US 



ft(eV) 



Q(eV) 



o.o 

1.51) 

1 75 
1 00 
0.75 
0.50 
0.25 
0.00 



FIG. 9. Calculated SXRS spectra from trans-NMA, in which both pulses are polarized parallel to the lab frame V axis. The 
pulses are Gaussian, 181 as FWHM in intensity, with center frequency set to either 401.7 eV (N) or 532.0 (O). From left to 

The two-color signals (bottom two rows) have their pulse 
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right we show the real part, imaginary part, and modulus of Eq. 
sequences given from left to right in chronological order, i.e. the ON signal results from having the O pulse come first and the 
N pulse come second. 



The contribution to the 2D signal coming from, for ex- ample, diagram 8; (see Eq. 13 1 is given by 



R (( a 2;gg") *«3;g"g' a l;g'g) fig 
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(<4 S + (r a , - ino 2 ) {w], g „ + (iv - 1^ 2 ) 2 ) 

) *a 3 .gU g ,a 1 . g ,g) (r g , (Ug,g + CJg,g„) ~ I (U g, g„ Q 1 + Ug'g^ 2 )) 
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Kg + (r fl - - »n0 2 ) K%„ + (r s , - in 2 ) 2 ) 



r 



and the other contributions follow similarly. As in the 
ID signal, we remove any zero frequency contributions 
(for fii and CI2) prior to plotting. 

Fig. [8] shows that the anti-Hermitian part of a, for 
both pulses is mostly diagonal, while the Hermitian part 
has significant off-diagonal contributions. Since the 2D 
signal depends only on the anti-Hermitian part of the po- 
larizability for the probe pulse, the probe acts as a filter 
through which the various valence-excited wavepackets 
created by the pump pulses interfere with each other. 



In Fig. [T0]the modulus of Eq. |22]is shown for the NNO 
pulse sequence, in which the first two pulses are resonant 
with nitrogen and the last with oxygen. Diagonal peaks 
are the largest features due to the large diagonal matrix 
elements of a 3 . The contributions Si and Sm contain 
diagonal peaks for g = g" . Likewise, in Sa and we 
get a diagonal peak when g 1 = g" . When different va- 
lence excited states are involved, off-diagonal peaks of 
two different varieties are found. Sa and Si v contain 
peaks at (^1,^2) = (w ge < , w ge » ) , i.e. u 2 is a valence ex- 
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FIG. 10. Simulated 2D-XRS spectrum from trans-NMA using 
an NNO pulse configuration, plotted as the modulus of the 
Fourier transform and separated into the contributions from 
the two types of diagrams in Fig. [2] The labels refer to the 
pulse center frequency and polarization of the three pulses 
ordered chronologically from left to right. In the NNO signal, 
the first and second pulses have their center frequency reso- 
nant with the nitrogen K-edge transition, and the third pulse 
is likewise tuned to the oxygen K-edge. Signals are plotted 
using an arcsinh nonlinear scale (see color bar) to highlight 
weak features. 



citation frequency of the system. In the other diagrams, 

fi 2 



ug'g' 



is the difference between system excitation 
frequencies. Inelastic peaks from the §j and Em contri- 
butions are therefore spectrally removed from the much 
stronger elastic peaks. Since each diagram contributing 



to the signal defined in Eq. 22 has inversion symmetry, 
further plots will omit the portion of the signal where 
fii < 0. 

Off-diagonal peaks are weaker than diagonal peaks be- 
cause the diagonal elements of ocj are much larger than 
the off-diagonal elements. To enhance these features, we 
plot the 2D signals using the nonlinear scale, 



arcsinh(S) = In (§, + \A + § 2 ) • 



(24) 



This scaling function interpolates between linear (for 
small §) and logarithmic (for large S) scaling and shows 
both weak and strong features (the modulus of the signal 
is taken prior to applying this scaling function). 

Plots of all 8 two-color 2D-SXRS signals are given in 



Figs. 11 and 12 Detailed interpretation of this informa- 
tion, relating the intensities of the peaks in a 2D spec- 
trum to the shape of the various orbitals and their un- 
derlying many-body wavefunctions will require further 
study. These signals may also prove to be a valuable 
test for the quality of electronic structure calculations 
on molecular systems. 

The ID spectra with the NO and ON pulse configura- 
tions are very similar, as may be seen in the bottom two 
panels of the right column of Fig. [9] However, the 2D sig- 



nals, especially in the spectral region 6eV > u>2 > — 6eV, 
show a much greater sensitivity to pulse combinations. 
Note the radically different spectrum that results in Fig. 
11 when the last two pulses are interchanged in going 
from the NNO to the NON configuration. 

In Fig. [13] we show the OOO spectrum from Fig. [12] 
on an expanded scale, together with several ID horizon- 
tal and diagonal traces. For comparison, we also show 
the corresponding traces from the OON signal using a 
blue dashed line. These two signals differ only in the 
frequency of the third pulse, the first and second pulses 
are the same. A trace along the diagonal line Cli = Q 2 
is shown in panel i and is very similar to the 00 1D- 
SXRS signal, the second row, right panel from Fig. [9] 
with narrower linewidths. Although the S2 (8.14 eV) 
is totally absent from the OOO signal, which is consis- 
tent with the ID results, this peak has amplitude in the 
OON signal. The matrix element of ao connecting So 
and S2, the square of which determines the amplitude of 
this peak in the 00 ID spectrum, may be small but it 
must be nonzero for this peak to show up in the OON 
signal. 

A horizontal trace, with a constant ^2 = 8.95 eV, indi- 
cating that the system is in a IS3) (5o| coherence during 
T2, is displayed in panel ii. The diagonal peak is domi- 
nant here, indicating a lack of off-diagonal peaks of the 
first kind, those arising from diagrams and The 
OOO and OON signals are largely the same, with the 
latter being larger in magnitude. Panel iii shows a sim- 
ilar horizontal trace, with Cl 2 = 8.14 eV corresponding 
to S 2 . This trace is much larger in magnitude when the 
probe pulse is resonant with nitrogen core transitions, 
consistent with the fact that the S 2 valence state is more 
strongly perturbed by a nitrogen than an oxygen core 
hole. 

The off-diagonal peaks seen in panels ii and iii result 
from the interference between diagrams 




where S n is the state giving rise to the flj frequency. 
The (fii,n 2 ) = (8.95 eV, 6.91 eV) peak in panel iii de- 
pends nonlinearly on the matrix element ao-SsS!- Only 
off-diagonal peaks carry information about the polar- 
izability between valence excited states. The diagonal 
peaks with contributions from all four diagrams in Fig. 
[2] are insensitive to these quantities. The diagonal peak 
at (Cli,Cl 2 ) = (8.95 eV, 8.95 eV) in panel ii, for example, 
depends on products like a 0l SoS3<xo;S3S3<xo;S 3 So where 
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FIG. 11. Simulated 2D-XRS spectra from trans-NMA, plotted as the modulus of the Fourier transform. The labels refer to 
the pulse center frequency and polarization of the three pulses ordered chronologically from left to right. In the NNO signal, 
the first and second pulses have their center frequency resonant with the nitrogen K-edge transition, and the third pulse is 
likewise tuned to the oxygen K-edge. Signals are plotted using an arcsinh nonlinear scale to highlight weak features. Traces of 
each signal along the diagonal are shown in red on top of each signal. 
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FIG. 12. Same as Fig. [TTJfor the other pulse configurations. 
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FIG. 13. (Left) An enlarged version of the OOO spectrum 
from Fig. [ll] plotted using a nonlinear scale shown on the 
color bar to the left. (Right) Horizontal and diagonal slices, 
plotted using a linear scale, of the 2D spectrum on the left 
(in red) plotted together with the corresponding traces from 
the corresponding OON (dashed, blue) to highlight the effect 
of changing the probe pulse in the three-pulse sequence. 



the polarizability between the ground and valence ex- 
cited states, which determine the linear SXRS signal, are 
multiplied by a diagonal element of the polarizability, 
seen to be much larger than off-diagonal elements in Fig. 
Ill 

Panels iv, v, and vi show diagonal traces where Q% 
is equal to the difference between two valence excitation 
frequencies. The peaks in panel iv result from interfer- 
ence of the following diagrams 



a3> 

«2> 




In panel iv the signal is much weaker with a nitrogen 
probe than with an oxygen probe though they are similar 
in shape. This is contrast with the two-color ID results 
from Fig. [T] where the Si peak is much stronger with a 
nitrogen probe than with an oxygen probe. 

Panels ii and v both show signals where the second 



pulse leaves the system in the S3 state, however they do 
so using the two different types of diagrams discussed 
above. The fact that these spectra are so different in ap- 
pearance demonstrates the wealth of information avail- 
able using 2D-SXRS. 



V. CONCLUDING REMARKS 

We have calculated and analyzed multidimensional 
stimulated x-ray Raman spectroscopy signals at the N 
and O K-edge. The signals are sensitive to the order in 
which the atom-specific x-ray excitations are used to per- 
turb and probe the valence electron configuration. The 
relative intensities of cross peaks could serve as a diagnos- 
tic tool for the comparison of different levels of electronic 
structure theory. By comparing frequency-domain RIXS 
and time-domain 1D-SXRS signals we show that the abil- 
ity to probe electronic dynamics at different atomic cen- 
ters is an advantage of the time-domain technique. 2D- 
SXRS offers a greater number of pulse combinations, and 
greater control over the preparation and measurement of 
electronic wavepackets. 

The 1D-SXRS signal can be decomposed into the over- 
lap of a doorway with a window wavepacket in Hilbert 
space. As discussed in Sec. |III| the 2D-SXRS experiment 
can be written as the sum of two valence wavepacket over- 
laps. More detailed analysis of this signal is needed to 
relate these terms to valence electronic dynamics. Stim- 
ulated Raman techniques in the x-ray frequency range 
will allow many new and interesting measurements to be 
performed. 2D-SXRS provides a new window into the 
correlated valence electronic dynamics made possible by 
new x-ray sources. 
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Appendix A: The Effective Transition Polarizability for 
Gaussian Pulses 



which can be simplified further to give Eq. [7] 



mi a i.- i u-i-i. u i r i. Appendix B: Rotationally averaged all-parallel 2D-SXRS 

The effective polarizability has a simple form when " v 3 B v 

the pulses have Gaussian envelopes. Assume the pulse s '§ na ' 
envelopes are 



Here we give the all-parallel 2D-SXRS signal as a ten- 

(Al) 

sor contraction over the polarization dependent expres- 
where aj is the temporal pulse width (equal to intenstity slon 

FWHM divided by 2-^/^2), which determines the band- ^zd/a j. \ ore u n i s /i ( iiD u nic u i 
width of the j pulse. Inserting this definition into the (Bl) 
time-domain expression for the polarizability (Eq. @, we where §2D ig thg sum of four terms (gee pig g 
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The Ti integral can be easily done by switching to polar 
coordinates, and invoking the identity 
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This integral can be written in terms of the complemen- 
tary error function, 
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with the 12th-order isotropic tensor^ 



Substituting Eq. |A6| into Eq. A5 gives 
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